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\ Abstract 

We consider a wide class of approximate models of evolution of singular distributions 
of vorticity in three dimensional incompressible fluids and we show that they have global 
smooth solutions. The proof exploits the existence of suitable Hamiltonian functions. 

The approximate models we analyze (essentially discrete and continuous vortex fil- 
aments and vortex loops) are related to some problem of classical physics concerning 
turbulence and also to the numerical approximation of flows with very high Reynolds 
number. Finally, we interpret our results as a basis to theoretical validation of numerical 
■ methods used in state-of-the-art computations of turbulent flows. 
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1 Introduction 

In this paper we consider the Euler equations describing the motion of an ideal (inviscid), 
incompressible and homogeneous fluid in M. n (n = 2, 3): 



^ + (u-V)u + Vp = 0, (1.1) 

and 

Vu = 0, (1.2) 

with the associated initial condition vi\t=o = uo(x), where u = (u , ...,u") is the velocity 
field, while the scalar p is the kinematic pressure. 

To solve (|1.1|) and (jl.2j) . the introduction of vorticity field is crucial. This quantity, that 
measures the "amount of rotation" inside the fluid, is responsible for the big differences 
between the 2D and the 3D case, see Wolibner (HE] and Yudovich 38 j. If n = 2, taking the 
curl of equation , one gets 

^ + (u-V)* = 0, 

where u := curl u = (0, 0, d\u 2 — c^u 1 ), that is a pure transport equation. Then, a bounded 
initial vorticity remains bounded for all the positive times. Let us denote by Xt(x) the 
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position at time t of the fluid particle that at time was in x. We have the following equation 
(equations for the characteristics) 

^^ = u(X 4 (x),t). (1.3) 

The vorticity field is simply transported along the trajectories of particles Xt, t > 0, remaining 
always perpendicular to the plane of motion. Beside, the vector u having two component 
identically zero behaves essentially scalar. 

When n = 3, the situation becomes much more complicated. The vorticity associated to 
u is a 3D vector field satisfying the following equation 

^ + (u • V) c3 = (w • V) u. 

In contrast to the 2D case, the vorticity u is no longer conserved and the fundamental 
complication arising in three dimensions is the presence on the right hand side of (u3 ■ V) u, 
the "vortex stretching term," which is -roughly speaking- of order of |cD| 2 . 

Some numerical simulations suggest that solutions of the ideal fluid equations may de- 
velop singularities and that the possible singularities may be detected by the blow-up of the 
vorticity field cD, see Beale, Kato, and Majda [2] and references therein. Moreover, accurate 
direct numerical simulations of 3D turbulent fluids show that the regions where |t3(x)| is 
big have an elongated "filament" shape, see for instance Bell and Markus 3 , Vincent and 
Meneguzzi [35] : and the reviews in Chorin [Hj and Frisch [17] . Thus, a consistent part of 
the current research activity is devoted to find suitable mathematical descriptions of these 
(and other geometric) vortex structures, see for instance Constantin In recent years also 
the use of the Lagrangian approach seems promising, to detect different intrinsic geometric 
properties of fluid flows, see Constantin |1U| . 

We briefly recall some recent contributions on the study of these filaments, motivated 
also by application to the numerical simulation: Chorin (§J proposed a model based on self- 
avoiding random walks, using ideas coming from the theory of polymers; Gallavotti |18j 
suggested the introduction of non-smooth curves (like Brownian motions), to avoid some of 
the divergences arising in the equations; Lions [2^j and Lions and Majda [26 proved a mean 
field result based on an approximate model of "nearly parallel" vortex filaments; Flandoli ^3] 
introduced a probabilistic description and overcame certain divergences in the energy by 
considering filaments with a fractal cross-section. Furthermore, by using this model, authors 
in JH] have defined certain Gibbs measures, while Bessaih and Flandoli 5 proved a mean 
field result. Recent overview on vortex computational method can be found in Cottet and 
Koumoutsakos [Ti] . 

On the other hand, relevant advances in applied mathematics have been obtained in the 
study of 2D models for fluids (like the quasi-geostrophic equations) that replicate the 3D vor- 
ticity behavior via "analytic and geometric analogies," see Constantin, Majda, and Tabak 
and Cordoba and Cordoba [12]. This is obtained with a smoothing or a modification (with 
appropriate scaling properties) of the Biot-Savart law. Finally, very recently Cordoba et al. 
(see ^3] an d references therein) studied the breakdown of smooth tubes of vorticity and other 
singularities concentrated on small but elongated sets. 

The basic goal of this paper is to show global existence of smooth solutions for models of 
line vortices, supplementing previous work in 0]. In particular, we shall exhibit a Hamiltonian 
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function. The main results will be derived by showing that the energy is conserved and also 
that the velocity induced by a singular distribution of vorticity is smooth enough. We shall 
consider rectifiable curves, but similar results can be also be proved for Holder curves and 
some rough paths. 

We shall also consider the vortex loop approximation (see Buttke [Jj and Chorin 8 ) for 
the equivalent equations of motion involving the impulse density. Also for this model we shall 
prove similar results of existence for large times of a smooth solutions. 

Finally, we discuss the evolution for approximate equations in which a smooth vortex 
filament is approximated by a discretization and we shall connect the existence of solutions 
with an appropriate choice or refinement of interior points. This latter model is motivated 
being the simpler numerical approximation to the line vortex evolution problem. 

Plan of the paper: The paper is organized as follows: In Section [21 we recall the motivation 
for the introduction of models with vorticity concentrated on singular sets together with 
classical methods to obtain integrable equations. Section introduces a proper Hamiltonian 
and gives the main analytical results: conservation of energy and estimates of the velocity in 
terms of the energy that prove global existence of strong solutions. Section uses the same 
technique to deal also with the problem of vortex loops in magnetization variables. Finally, 
in Section a question connected with the resolution needed to simulate the behavior of a 
line vortex is put in the framework of the previous results. 

2 Vortex motion and approximate models 

To introduce the problem we shortly recall the discussion in Section 2 of [I] , where the reader 
can find full details. Curves that are always parallel to the vorticity vector are known as 
vortex lines, and vortex lines passing through the points of a closed curve define a volume 
called vortex tube. 

A vortex filament is a vortex tube that is immediately surrounded by irrotational fluid, 
even if the usage of this expression is neither uniform nor consistent. It can be also used to 
design infinitesimal vortex tubes or also patterns of vorticity that are not properly formed by 
vortex lines. 

A line vortex is a singular distribution in which infinite vorticity is concentrated on a 
line, such that the circulation around a closed circuit threaded on by the line is finite. A 
line vortex is the result of a limiting process in which a vortex filament (of finite strength) is 
contracted to a curve, the strength T > being kept constant. 

The laws of vortex motion, first formulated by Helmholtz [20], state that in an ideal fluid 
vortex lines and tubes move with the fluid and that the strength of each vortex tube (flux 
through any cross-section) remains constant in time. 

Hence, we consider an ideal fluid with the vorticity u3 concentrated on a smooth curve 7. 
Neglecting the transverse size of the filament, we can write the vorticity field u generated 
by 7 (that is described parametrically as a continuous function 7 : [0, 1] — > K 3 such that 
7(0) = 7(1)) and, formally, 

u;(x, t) = T f <5(x - 7(£, t))i^, t) d£ Vx G M 3 , Vt £]0, T[, (2.1) 
J 

where <5( . ) is the Dirac's delta function, £ is the arc-length, the parameter t represents the 
time, and 7g is the derivative with respect to the arc-length. 
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A discussion on the precise meaning of the formal vorticity field (|2.1I) and its distributional 
rigorous definition, can be found, for example, in the appendix of 

Remark 1 The above formula can also be written formally as the distribution 



w(x,i) = T / <5(x-7 5 )^7, 
Jo 



meaning that the integral could be the Riemann-Stieltjes integral or also that included in a 
more general setting. In the sequel we shall not write explicitly the dependence of u> and u 
on 7 since we shall always consider vorticity (and corresponding velocity) "generated" by a 
distribution concentrated on a curve 7. 

It is well-known that the kinetic energy associated to (|2,lj) (in the sequel we shall denote by 
(a, b) the scalar product in M 3 of the vectors a and b) 



E{t) 



r 



8vr7 Jo l7(£,*) 



7(77 t)\ *)'%(^' *)) d Z dr l> 



is infinite for any reasonably smooth curve 7. This explain the reason for the study of 
approximate models (as the Localized Induction Approximation (LIA) or the Rosenhead ap- 
proximation, see Saffman [3S]) or for the same problem with curves supported by non-smooth 
lines (see Gallavotti |18j and recent results in 6 ). To derive the evolution equation satis- 
fied by a line vortex the essential tool is the Biot-Savart formula that gives a representation 
formula for the velocity, in terms of the vorticity: 

u(x) = --3- f (J^ A d5(y)dy, 
4vr J R3 |x - y| J 

where a A b € K 3 is the vector (exterior) product of the vectors a, b € M. 3 . Inserting the 
expression (|2.1[) of the vorticity field in the Biot-Savart formula, using the equation for the 
characteristics ()l,3j) . and imposing that the curve 7 is transported by the velocity field u, we 
get the following equation 

This equation (relation between its solutions and weak solutions to the 3D Euler equations can 
be found -for instance- in the appendix of 0J) has been used for the first time by J.J. Thomson, 
who proposed a cut-off in the integral in order to control an infinite energy, see |33j . An 
approximation based on the Taylor expansion of the kernel in the singular integral gives the 
Localized Induction Approximation (LIA) solved by Hasimoto while recent results in 
this direction are those by Klein and Majda |23l I24| . 

The de-singularized model studied analytically in |1] is that proposed by Rosenhead 32 , 

07 _ r r 7K,*)-7M) 



dt 4»i [(-J( f ,t)-7(,,t)) 2 + M 2] 3/2 

for some [i > 0. The choice was motivated by the fact that (|2.2j) gives advantages for numerical 
integration and, starting by Moore [HO], it has been used in the applied context for some 
calculation related to aircraft trailing vortices. 
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2.1 On the Rosenhead model 

Regarding the Rosenhead model, a first analytical result has been recently proved in 0j. In 
that reference is shown (as a particular case) that if the initial curve is smooth and closed then 
the initial value problem is well-posed in natural Sobolev spaces. In the sequel we shall use 
the customary W k,p and H k = W k ' 2 spaces, see Adams pQ. To our knowledge the following 
is essentially the first rigorous existence and uniqueness result, known for this model. 

Theorem 1 [See ^] Let 7 (£) £ fli(0, 1), where iJi(0, 1) denotes the subset of closed 
curves R h>1 3 belonging to the Sobolev space i? 1 (0, 1). Then, there exist a strictly positive 
T* < T and a unique curve 7 such that 

a) 7 € W^iO, T; L 2 (0, 1)) D L°°(0, T; H^(0, 1)); 

b) 7«,t) =7o(0-^ f f „ ; ^' T) ~^ r) A7,(r)^r, 

4Wo 7o [(^(e,r)-^,r)) 2 + M 2 ] 3/2 

£/iat is a strong solution to (j2,2j) in the time interval [0, T*). 

In addition the time T* > may 6e estimated from below and it depends just on the 
H 1 -norm of the initial datum. 

A similar result has been recently proved in a more general stochastic context in [g] with a 
different target and weaker hypotheses on the smoothness of the curve 7. Anyway, in both 
papers the question of the global solvability left open and some continuation criteria involving 
the length of the curve itself were derived. 

Our main goal is now to show that, for a class of regularizations (including the one above) 
the above local theorem is in fact global. 

2.2 A class of regularized problems 

An interesting class of regularized evolution equations can be defined as follows. Let <p : M 3 — > 
R be a scalar function and define the velocity field u associated to the vortex filament 7 as 

u(x):=(curl^)*cD= / V^(x - 7(e)) A <^(£) d£, (2.3) 

J 

where the convolution * has to be intended in the sense of distributions. Formally, this 
formula is such that curlu = u3. 

On the kernel <p we impose the following conditions that will allow the energy to be 
well-defined. 

Hypothesis A. In the sequel we shall assume that: 

- The function ip is even: 

¥>(-*) = p(x); (A.l) 

- The function tp has a real and non-negative Fourier transform: 

£(k):=/ e i{k '^(p(5c)dx> VkeR 3 ; (A.2) 
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The Fourier transform £>(k) is integrable over R 3 : 

i dk < +00. (A.3) 

- The function ip is smooth enough, in order that 

/ (1 + |k| 2 ) V(k) dk < 00; (A.4) 

These assumptions are satisfied by a wide range of function 92, covering also well-known cases 
with a physical meaningful interpretation. 

Remark 2 (An explicit relevant example) The function 

r 



IxP + m 2 ) 1 / 2 



for some \i > would be a natural candidate satisfying all the above assumptions. Assumption 
(A.l) is trivially satisfied. Regarding the assumption (A.3), the Fourier Transform may be 
expressed in terms of a Bessel function of the second kind y(z) = Bessel(l, 2), defined as 
solution of the differential equation z 2 y" + zy' — (z 2 + l)y = 0. The condition on the non- 
negativity of the Fourier transform can be verified by a direct computation: 

r e »(k,x) 

r*t e " R ^;f P "" J <* 







00 p -\k\y(4tyt^ 



where we used the fact that 



dt > 0, 



A -l/2 =7r -l/2 / e -Xt t -l/2 d ^ 

Jo 



and employed the explicit expression of the Fourier transform of the Gaussian kernel. 

Regarding (A.3) the integral over M 3 can be split into the inner and outer part. The first 
one is bounded since the function ^i?(k) is 0(|k|~ 2 ) near the origin, hence the inner integral 
converges. On the other hand, the decay at infinity necessary to show the convergence of the 
outer integral comes directly by observing that ^i?(k) = 2fi 2 Bessel(l, 2|k|/i)( / u|k|)~ 1 . The fast 
decay at infinity (enough to have both (A.3) and (A.4) satisfied) derives from the properties 
of this special-function. 

In this case we obtain, by using as kernel the function ipn, exactly the Rosenhead model 
where the velocity field is given by the equation: 



fi(x) = r ! l ; A 7^(0 d£- 

Jo (|x-7(OI +^ ) ' 
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The evolution problem can be set up in different Banach spaces of closed paths where local 
solution exists and are unique. For example we can have solution for initial condition in the 
Sobolev space iTi(0, 1), in the space of Holder continuous with exponent greater than 1/2, 
and also in some spaces of rough-paths, see jUIH]. 

Details of the proofs (and some specific assumptions) depend however on the functional 
setting, so we will show the proofs in the case of solutions living in H^(0, 1), using the 
framework of Theorem 



3 Energy and global solutions 

For the evolution problem 

§p&t)=um,t),t)- (3-i) 

with velocity given by ()2.3|) we can identify a function which plays the role of a "kinetic" 
energy. 

Definition 1 (Kinetic energy.) The function below Tt : y i— > R, defined on the space of 
smooth curves, is the "kinetic energy" for the smoothed evolution problem associated to (|3.1|) 



1 



l ri 
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n ^ : =2l- I- ^® ~ 7fa))<7e(0, 7„fa)> d^drj, (3.2) 



The main result we shall prove is that H(y) is constant in time. Then, the existence of the 
energy allows to exploit a-priori estimates to have global existence for initial data with finite 
energy. 

Remark 3 This particular expression for the energy has been discussed also by Marsden and 
Weinstein and Holm \21^ . they show that on the space of closed curves in R 3 there exists 
a natural Poisson structure and that, formally, with respect to this structure the function h'J.ty) 
is the Hamiltonian function which generates the flow described by I.V.i)) . 

First, we show that the function 1~L(ff) is well-defined. 

Lemma 1 For each smooth enough curve y, it holds that < H(y) < +oo. 
Proof. In terms of Fourier variables, the energy 7i(ff) can be written as 

H(y) = ^3 £ £ ^(k)e-**-7W) i^)%(v)^d V 



1 



2(27r) J 



£(k) 



*<k,7(e)> 



dk 



which, thanks to the assumption on the non-negativity of the Fourier transform of ip, proves 
that the energy is a non- negative quantity. Moreover, for any y E Hi(0, 1) we have the 
obvious estimate ^ 

l7 e (0l^, 



e*<k,7(0>^(£)^ 



< 
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actually we observe that sharper estimates are possible, even if they are not needed at this 
stage. Finally, we have proved that 



The final observation is that both integrals are finite (use assumption (A. 3) and the fact that 
7 is rectifiable a curve), showing that ^(7) is a well-defined energy for any 7 € H^(0,1). □ 

Next, we show that the function behaves as a constant of motion if 7 evolves under 

the flow associated to u. 

Lemma 2 Let~f{t,£) be a local smooth (as those of Th.^\) solution of the problem then 



Proof. We shall show by an explicit computation that the energy ^(7) is invariant and 
the completely anti-symmetric Levi-Civita tensor will be used to write, with the Einstein 
repeated indices convention, that 




dnm,o) 



= 0. 



dt 




Then, with explicit vector notations, we have: 



dnm,o) 



1 



L 



1 



dt 



+ 




where, for simplicity, dtjit) := d~f{t)/dt. 

By integrating by parts the ^-integral in the second term we get 




where we used the fact that 




Next, by definition of u we have 




= 
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where we used in sequence the equation (|3.1|) of motion for 7(f) and the complete anti- 
symmetry of the tensor e^. □ 

The next step is to show that if the kinetic energy is bounded, then the velocity u associ- 
ated to the evolution problem is smooth. In particular, it will follow that the velocity induced 
by vorticity concentrated over a i? 1 -curve is very regular. The regularity of the smoothing 
kernel is inherited by the velocity u even if the framework is that of a singular problem. 



Lemma 3 For any < n £ N, we have the bound 

1 



V n u 



< 



2tt 3 / 2 



k| 2 ( 1+n ^(k) dk 



nl/2 



^ 1/2 (7), 



(3.3) 



provided that the integral J K3 |k| 2 ^ 1+n ^(k) dk is finite. 

Proof. The proof follows easily by using Cauchy-Schwartz inequality. In fact, it follows that 



1 



u 



< 



(2tt) 3 
1 

(2^)3 

1 

2vr 3 /2 



e*<M>£(£)(i£) A / e-^^HOd^dk 
Jo 

1/2 



\k\ z p(k) dk 



|k| 2 £(k) dk 



1/2 



£(k) 
^ 1/2 (7), 



7(0 dti 



lV2 



dk 



and in a similar fashion we can bound all the derivatives of u, provided that ip is smooth 
enough: 



|V n u| 



< 



1 



2tt 3 / 2 



|k| 2 ( 1+ ">^(k) dk 



1/2 



^ 1/2 (7)- 



□ 



Essentially we shall use this lemma just for n = 1 to prove global existence of solutions, 
hence assumption (A. 4) will be enough. On the other hand, to prove higher regularity of u 
these extra conditions on tp are needed. 

Theorem 2 The evolution problem (|3.1|) has a unique global solution for any initial condition 
in H^(0,1). 

Proof. Denote tp G H\{Q, 1) the initial condition, then we already know that local solution 

7(i,£) exists in a (possibly small) time interval t S [0,T*[, depending on the L 2 -norm of ip^, 
recall Theorem ^ However, by taking the derivative of (|3.1|) with respect to £, and with 
and integration over [0, t], we have now the additional a-priori estimate (|| . || denotes the 
L 2 (0, l)-norm) 

||7{(*,0ll<ll7{(0)||+ / ||Vu( S )||||7 5 ( S ,6ll^ 

Jo 



\7t(o)\\ + CHl/2 tf(°)) WM s '0\\ds, 
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where the fist line is an easy bound for the norm of the solution, the second line comes from 
Lemma|21 and the last line from the fact that the energy is constant along solutions (as stated 
in Lemma |2I) 

Then, by using the Gronwall inequality we have 

sup ||7 e (t,6ll<ll7 5 (0)||e cwl/2 ^°)) T , 
te[o,T] 

which guarantees existence of global solution, since any local solution, having the i? 1 (0, 1) 
norm uniformly bounded can be uniquely continued up to any positive time T. □ 

3.1 Holder initial conditions 

In |2j the authors prove that for sufficiently regular functions ip the vortex filament equa- 
tion has local solutions for initial conditions 7(0) which are Holder continuous functions of 
exponent a with a > 1/2. In this case, the line integrals of the form 

(f(7(t,0), 7{(U)K = /V(7(*,0), da(t,0) (3.4) 

J o 

must be understood as limits of Riemann sums: such limit exists if the function f : IR 3 — > M 3 
is at least C 1 and this process defines the integral in the sense of Young |37| . 

Moreover, in the same paper, the authors extend the local existence result, to a class of 
initial conditions living in a space of rough paths (see the book of Lyons for more details |271 
28 ) that are special class of Holder paths with some additional structure. For rough paths it is 
possible to define the line integrals (as those in eq. (|H.4|) ) by means of natural "renormalized" 
Riemann sums. In this way it has been proved that the filament equation (with sufficient 
regularity of the kernel tp) admits local solutions starting from almost every 3D Brownian loop 
(i.e. curves chosen accordingly to the Wiener measure restricted to the subspace of closed 
curves) . 

The proof of the constancy of the energy extends also to these different functional setting 
and the existence of global solutions can be proved along the same strategy used for H 1 paths. 
We leave the technical details to the interested reader, since they are outside the scope of the 
present paper. For example, the extension of Lemma[2]to Holder solutions can be done using 
convergence of approximating discretizations and the computations contained in Sec. [5J 

Remark 4 The relevance of this additional result can be seen in the light of the K41 theory 
and the Onsager conjecture on possible singularities in the sense of C°' a velocity fields. This 
may suggest that velocity behaves as Holder continuous velocity field (see \lTJj ), hence as 
a singularity diffused all over the fluid. More recent developments link intermittency with 
concentration of singularities on small sets, described by Holder -hence (multi-) fractals- sets, 
see {T^ E3/. 

4 Impulse formulation 

In this section we consider another singular distribution of vorticity that is well-known and 
useful in the scientific computing, since it involves vector fields (as the "magnetization" ) that 
are possibly with non-vanishing divergence. In particular, we shall consider Buttke loops that 
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(according to the terminology of Chorin jHJ) are small loops of vorticity in an irrotational 
background, which evolve according to 3D Euler equation. The first study of the kinematic 
interaction of an immersed body in an inviscid irrotational flow dates back to Kelvin |22] and 
his analysis was based on the study of fluid impulse. It was in that paper that he introduced a 
model of "core-less vortices" in order to explain some experimental facts on eddy formation at 
the boundary. The interaction of a vortex ring (line) in an inviscid flow has been also studied 
by Roberts [HJ) by using as canonical variables the position of the centroid of each ring and 
its impulse (called by Roberts "momentum of vorticity"). A proper Hamiltonian formulation 
of the 3D Euler equations have been discovered independently by different authors, see for 
instance Osedelets |34j for a development in a continuum setting, with canonical variables 
the position and the impulse density. Then, around 1990 Buttke 7^ linked the discrete 
formulation of Roberts to fast and efficient discrete methods for the numerical simulation of 
turbulent flows. We consider also this setting, since by using essentially the same techniques 
of the previous section we are able to prove a global existence result also for the latter model. 

4.1 Buttke loops and a discrete problem 

The main idea (refer for instance to |H]) is to introduce a new variable m (that is called 
magnetization or vortex magnetization), that is obtained by adding to the velocity u a gradient 
at t = 0: 

m = u + Vq. (4.1) 

The unknown m does not satisfy the incompressibility constraint, but it is with compact 
support and it satisfies 

curl u = curl m. 

Then, we can interpret rri as essentially local, while Vq is an extensive field. The decompo- 
sition H4.1|) resembles the Helmholtz decomposition, even if m is different from u, or from 
momentum density. The vector m is related to the so-called "effective vorticity" (compactly 
supported) by the relation £ = curl rri. Finally, the introduction of the new variable m can 
be seen as a gauge transformation, known as "geometric gauge." 
The equation of evolution for rri can be easily determined: 

^ + (u ■ V) m + m • (Vu) T = 0, 

where u = Pm, P denoting the Leray projection over divergence-free vector fields. The main 
point is that we have an equivalent equation for u plus an arbitrary gradient at the initial 
time. A proper choice of this gradient leads to the study of vortex loops. In fact, if curl u has 
support within a small ball B, then the resulting m has support within the same ball and 

m(x) = M p(x — xg) for some xg £ B. 

Here M is a vector in M 3 , while p is a non- negative smooth function with support within 
B and such that J K3 p(5c) dx. = 1. The "magnet" Mp(x — xg) has a simple interpretation: 
the velocity field induced is the same of a small vorticity loop (or vortex ring), with M 
perpendicular to the plane of the loop and |M| = TnR 2 . Then, from the physical point of 
view the vector m can be interpreted as a vortex dipole density. 

The above interpretation can be used to deduce a Hamiltonian system for a more com- 
plicated distribution of vorticity, see Buttke 7 . This approach works as a de-singularization 
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of the problem: the resulting discrete dynamic we shall consider includes in an implicit way 
a regularization that will allows us to prove global existence of smooth solutions as well as 
good computational properties investigated in jjj. 

The discrete model is derived by considering a finite sum of magnets and denoting by x ft , 
for a = 1, . . . , N, the positions of the loops and by m Q be the corresponding magnetization 
vector, that is not anymore a constant vector. Then, the velocity field u of the fluid is given 
by the solution of the equation 

u = M + Vq, 



with 



N 



M(x, t) ™a(t)p(x. ~ X c 

a=l 

where the scalar <p is determined by the equation 

Aq = -V-M. 
The equation of motion for M is (component-wise) 
dM* 



dt 



+ ui V i M i = - MP VV 



i = 1, . . . ,3 



and the "loop particles" moves according to the system of ordinary differential equations 

3 

W.VVfxJ 

(4.2) 



m„ = -^m£vV'(x c 

3=1 



x„ = u x f 



where the index a run over particle labels. 

Equations (|4.2j) form a Hamiltonian system, if we consider the conjugate pairs of variables 
(xq,, m tt ) a= i jv with symplectic structure 



N 3 



a=l i=l 



and Hamiltonian function 



1 N 

H({x a , &a}) = ^ m Q • p(x Q - x^) + (m a • V)(rii /3 • V)$(x Q - x^), 
a,P=l 

where <3? is the solution of A<I> = p. Hence system (|4.2|) takes the canonical form 

dH 



x n 



dH 
<9x„ 
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Remark 5 The analogy between magneto-static and fluid mechanics is that to the electric 
current it corresponds the vorticity and to the magnetic induction it corresponds the velocity. 
Since the name "magnetization" may be misleading, we prefer to use Chorin's notation and 
to refer to them as Buttke loops or vorticity loops. 

In order to show that the energy is constant on solutions we pass to the wave-numbers 
notation. Let p(k) be the Fourier transform of p(x). The energy H({x a , m a }) can be written, 
in terms of wave- numbers, as follows 



H({x a ,m a }) 



1 



2(2vr) J 



p(k) 



N 



i(kx Q > 



a=l 



dk, 



where II - is the projection in the plane orthogonal to the vector k. From this representation 
it is clear that the energy is positive definite. Moreover, the velocity field has the form 



1 f N 



i<k,x-x„) d £ 



and the tensor Vu appearing in the equation for m is 

1 f N 
[Vuf 3 (x) = I p(k)(-ik) 1 Y,( U C m «) j e~ i{k ' x - Xa) dk. 



(2vr) J 



a=l 



4.2 The periodic problem 

In order to avoid problems related with the decay at infinity and possibly non-converging 
integrals, we restrict to the periodic setting. Thus, we fix a box of linear size 2tt and consider 
the periodic version of the above problem, i.e., we substitute Fourier series expansion to the 
Fourier transform. 

The energy is now given by 

2 



W({X m a }) 



1 



2(2tt) 3 i 

kez 3 



N 



m Q e 



i(k,x«) 



a=l 



where k runs over Z 3 and the coordinate x a is restricted to the box Aj, =] — vr,7r[ 3 . The 
velocity is then 



u(x) = 



(2vr) 3 i 



N 

E p& n k E^' 

ke» a=1 



-i(k,x-x a ) VxG]-vr,7r[ 3 . 



By using standard techniques we are able to prove the following result. 

Proposition 1 The solution for the evolution problem of a finite number of loops in the 
periodic box A l exists ( and is unique ) for any positive time, provided that the initial condition 
has finite energy and the function p(k) is non-negative and satisfies the decay estimate 

|k| 4 p(k) < +oo. 

kez 3 
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Remark 6 The condition of decay for |k| — > +00 required in Proposition^ is easily satisfied 
if p(x) is a smooth enough function. The condition on the non-negativity of jo(k) is not 
automatically satisfied for each smooth function p(x). A wide class of smooth functions with 
compact support and non-negative (hence real) Fourier transform can be identified as follows. 
Consider a smooth (say C°° ) "bump" function b(x) over the ball B(0,5) and, in addition, 
suppose that b(x) is "even, " implying that the Fourier transform has vanishing imaginary 
part. Then, define p(x) := (b* b)(x) and this will turn out to be smooth, non-negative, and 
null outside the ball B(0,25) of radius twice that of the original one. Finally, by using the 
convolution theorem it will follow that 

p (k) = bTb (k) = b(k) • b(k) > 0. 

Proof of Proposition ^ Again the proof is based on an energy estimate and an a-priori 
bound of the velocity in terms of the energy. The main point of the proof is the following 
Cauchy-Schwartz inequality: 



ax < 



1 



(2vr) 3 i 

kez 3 



N 



m„ e 



i(k, x-x a ) 



a=l 



< 



1 



E 

kez 3 



1/2 



kezs 



N 



m a e 



i(k, x-x Q > 



0=1 



1/2 



2tt 3 / 2 

and in the same way it easily follows that 



E ?S) 

kez 3 



1/2 



H 1/2 ({x a ,m a }). 



lull + llVull + ||V 2 u|| < 



1 



2vr 3 / 2 



E IklV(k) 

kgz 3 



1/2 



H 1 / 2 ({x Q ,m a }), 



showing that we can control the regularity of u in terms of the energy Tl({x a , m Q }). By 
a direct differentiation it is also clear that 7Y({xq,, m Q }) is constant along (smooth enough) 
solutions of the equations Q4.2j) so that we have good control over u for any time. This implies 
the global existence of solutions. □ 



5 A class of discrete models of evolution of vorticity 

In this section we study a class of discrete models for the evolution of vorticity which includes 
discrete versions of the line vortex model. These models could be of interest in the numerical 
simulation using lattice methods as suggested in [S] . In the sequel we shall always assume on 
the function (p the four natural hypotheses (A.1)-(A.4), we introduced in Section f2. 21 
We consider distribution of vorticity giving rise to a velocity field of the form 

N 

u(x) = E V ^ - *«) A 

a=l 
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that is parametrized by the set of 6N variables {Xq,, ^ a } a =i,...,N- The physical interpretation 
of this equation is that of describing the superposition of N vortex "blobs" (i.e. concentrations 
of vorticity) indexed by a, situated at points x a , and whose vorticity is directed along the 
vector £ a . In this case these blobs determine the initial datum and there are no assumptions 
that the solution persists being of the same type for all positive times. This is the difference 
between blobs and loops considered in the previous section. 

A priori we can think to each of the variables {x a ,^ a } a as an independent degree of 
freedom. However, if we set £ a = x Q+ i — x a (with the understanding that a + 1 = if 
a = N), then this model corresponds to a natural discretization of the line vortex model 
previously described. 

For the moment we do not impose any constraint on the variables £ (either static or 
dynamic) and we just impose that the points x Q are transported by the flow u so that we 
have 

x a = u(x Q ), a = l,...,N. (5-1) 

This is the discrete counterpart of (|1.3|) . If we want to address the problem of the global (in 
time) dynamics of these models we can introduce the following useful quantities 

N N 

C = Y,\L\ and A = Y,\L\ 2 

a=l a=l 

that are respectively the length and the quadratic variation of the curve itself. Moreover, in 
analogy with the Hamiltonian introduced in Section |21 we define an energy function (which 
controls the magnitude and regularity of the velocity field u) as follows: 



1 N - - 

ft({x Q ,£ Q }) = - E V(*a-*0)(5oi«/j}- 



a,0=l 

It is not difficult to show that we have an analogous of Lemma EJ which allows us to 
control u in term of 7i and, of course, we have the following trivial bounds 

||V n u||L°° < c n>ip C, and Ti < ||v?||z,oo£ 2 , 

where the constants c n ^ depends only on ip and the number n of derivative (we assume that 
(p is sufficiently regular for all the constants in this section to be well defined. For example it 
is enough to take <p infinitely differentiable.) 

Let us now compute the time derivative of the energy. Taking into account eq. (j5. 1|) . but 
leaving the chance of an arbitrary time-derivative dt£ a for £ Q , we obtain 

dH{{ ^ a}) = E E a^-v^-x^)^,^). 

' a,/3=l a,/3=l 



We now analyze the second sum appearing in the above equality, in order to recover some 
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symmetries or cancellation properties. We get 

N N 
1 = E ' V <p(*a - X»(£ a , ip) - Co, • u(x/9) C/3 " V<p(x a - x» 

a,/3=l a,/9=l 
N 

= E A v ^(*° ~ §3 A "^)) 

a,/3=l 
AT 

0=1 

and this implies that we can write the following equivalent expression for the time derivative 
ofW({x«,£ a }): 



dH({x a ^ Q }) 
dt 



N 



N 



Consider the second term and write 

N 



a,/3=l 
A? 



Af 3 



/3=1 a,/3=l i j'=l 



where 



$(x) = ^(x-x Q )(£ a , u(x)>. 



o=l 

Finally, the rate of change of the energy takes the form 

N 3 



d_ 

dt 



N 



H({* a ,L})= E E^-^)^[^- vV (^)^]+E^-v$(x^). 

0,0=1 i,j=l 0=1 



We can rewrite the above expression by setting \&(x) = Y2 a =i ~ hence obtaining 



d_ 

dt 



3 3 



N 



W({x Q ,U) = E [ d ^ - ^W^] +Y,£p ■ V<&(x», (5.2) 

0=li,i=l ' 0=1 



with = curl\I/ since, by definition, u(x) = curl^E'(x). By using a first order Taylor 
expansion with integral remainder, the introduction of the auxiliary function <3? allows us to 
rewrite the second term of (|5.2|) as follows: 



N 



N 



N 



(3=1 P=l f3=l 



(5.3) 
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where is the matrix 

Jo 

We observe that analogously to Lemma|31we can control the function * in terms of the energy 
7i according to the following lemma, whose straightforward proof is left to the reader. 

Lemma 4 For any < n £ N, we have the bound 

In the same way we can bound <I> since, e.g., it holds 

||$||l°° = ||* • curl#|| L oo < Cip H, 

where is some positive constant depending just on the Sobolev regularity of ip. 
At this point (see also [S]) we may recognize to have two interesting cases: 

A) The discrete filament model, obtained by setting £ Q = x a+ i — x a . In this case the first 
term in the r.h.s. of Eq. ()5.3j) vanishes due to the fact that 

N N 

+ Ip) - $(x»] = J>(x> +1 ) - <D(x»] = 0. 

(3=1 13=1 

B) A blob model, obtained by considering each couple (x a ,£ Q ) as describing (an approxi- 
mation of) a vortex blob and for which the dynamics of the £ a is fixed by prescribing 
that these vectors are transported by the flow at the point x a , i.e. 

i a = (Vu(x a )f.t 
In this case the first term in the r.h.s. of Eq. (|5.2|) vanishes identically. 
For these two models we prove the following result: 
Lemma 5 In the discrete filament model we have 

\d t H\ < Cl HA + c 2 H 1/2 A. 

In the blob model we have 

\d t H\ < c 3 NH 1/2 + c 4 H 1/2 A \d t H\ < c 5 £H 1/2 . (5.4) 
Moreover, in both models we have 

\d t A\ < c 6 H 1/2 A \d t C\ < c 7 7i 1/2 C. 
The positive constants Ci, for i = 1, . . . ,7 , depend only on the function ip. 



|V n *| 



< 



kj 2n £(k) dk 
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Proof. Given Lemma^J the only nontrivial part is to justify the first term appearing in the 
bound for the discrete filament model, i.e., c\HA. This term comes form the expression 



JV 3 

p=i i,j=i 

appearing as the first term in the r.h.s. of eq. (|5.2|) . In the discrete filament model we have 
(note that £ depends just on t, hence the partial derivative is in fact a total derivative) 

£ a = x a+ i - x a = u(x a+ i) - u(x Q ) 
and, by a first order Taylor expansion, we get 

L = [Vu(x Q )] T • £ a + e a (£ a ® u, 

where the tensor Q a is given by the formula 

= [\i-p)v^a k (5i a + pi a )dp. 

Jo 

Then, we can bound G as follows: 

JV 3 JV 

\G\ < \\*\\ L ~ £ Yl l 5 4 - W(x» 4l = 2 i^fe ® 4)1 

/3=1 i,i=l J3=l 
JV JV 

< J^il^psupie^'*! < ii^ii^iivuHloc^i^i 2 



/3=1 



J3=l 



The other bounds are proved similarly. 



□ 



Remark 7 PFe noie i/iai the term c^NTi 1 / 2 in the rate of energy change of the blobs model 
derives from the trivial bound 



N 



13=1 



< 2jV||$|| L o 



which, however, misses any possible cancellation property of the various terms. Heuristically 
we expect that, if the evolution is constrained in a bounded volume (e.g. by imposing periodic 
boundary conditions on a torus), the distribution of the set of points {x a } a should not differ 
significantly from that of the points {x Q + £ a } a and thus that the two sums ^ a $(x a ) and 
$(x„ + £ a ) should not be too different (at least not of the order of N). 
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5.1 Heuristics about adaptive numerical algorithms for vortex filaments 

The results of the previous section suggest a strategy to implement a numerical algorithm 
which approximates the global solution of the filament equation. Here we do not pretend 
mathematical rigor but only lay down some possible developments of the main results of this 
paper. 

So let us imagine to start a numerical approximation of the vortex model with some 
arbitrary discretization level (with iV points) and follow the solution up to the instant when 
the variable A(t) is larger than a fixed value, e.g. 1. At this point proceed to a refinement 
of the discretization up to level M > N in order to ensure that A(t) become lower than 1/2 
and continue the simulation. Upon increased discretization the energy change, but it is not 
difficult to argue that this change will be small (since in the limit of infinite discretization 
the energy will converge to a limit). So in first approximation we will neglect this change 
(actually is possible to devise a discrete model where the filament is replaced by a piece-wise 
approximation, in which, upon refinement, the energy does not change at all). In this way 
we ensure that, at any time A{t) < 1 and so we have also an (approximate) bound for the 
growth of 7i: 

H(t) < n(0)e m!lK{ci ' C2}t 

which implies that there will be no blow-up of the solution. Of course since \dtC\ < cTi^^C 
there cannot be blow-up for the length of the discrete filament as long as the energy remains 
finite. 

This algorithm is very similar to the heuristic techniques that are already adopted in nu- 
merical simulations of vortex filament dynamics, but here we can identify a clear quantitative 
criterion when to perform the refinement. 

Recall also that the rate of change of energy is, in this case, given by 

, N N 

-w({x a ij) = ^(#(x /3 ),G /3 ^0^) + ^(^,e^ /3 ), 

(3=0 (3=0 

cfr. eq. 1)5. 2|) and Lemma El Then, a more refined algorithm to decide where to increase the 
resolution of the discretization is to find points along the filament where the quantity 

j p = (*(x», ® Zp) + dp, Qpip) 

is large. These points are responsible for large positive contributions to the change of energy 
and thus can be associated to points of instability of the finite-element approximation of the 
dynamics. By refining the mesh only in the neighborhood of these points we diminish the 
rate of change of the energy and control the discrete evolution. Currently we do not know if 
the latter refinement technique is too expensive from the computational point of view. 

Finally, we observe that the presence of N in the first term on the r.h.s of (|5,4j) prevents 
to obtain the same result for the vortex blob model, since in this case refinement will enlarge 
this term. 

6 Conclusions 

We considered several classes of problems related to the evolution of special patterns of 
singular vorticity. Various kind of regularization have been analyzed and they are based 
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either 1) on a smoother version of the Biot-Savart law, linking velocity and vorticity or either 
2) a "point vortex" special expression of the solution. Under reasonable assumptions we have 
been able to prove global-in-time existence of solutions for these various models, employing 
the Hamiltonian structure of the underlying 3D Euler equations. 

In the case of the vortex filament we have also been able to consider curves that are 
not very smooth, solving an approximation of the line vortex equation containing also the 
classical Rosenhead approximation for this problem. In the case of small vortex loops we 
solved a periodic problem involving a finite number of them. 

In the final section we used these results to obtain, as by-product, some results on a 
couple of discrete models. The discrete models we have been able to analyze are a discrete 
filament model and a model of 3d vorticity blobs. For the latter we have been unable to find 
global existence results, while for the former we linked global regularity of the limit problem 
(smooth curve) with a numerical procedure to properly simulate this vorticity pattern. In 
particular, a quantitative refinement technique has been suggested, in order to produce an 
effective numerical simulation for arbitrarily large times. 
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